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We show that there exists only a quite limited number of higher dimensional operators which can naturally 
lead to a slow decay of dark matter particles into monochromatic photons. As each of these operators inevitably 
induces decays into particles other than photons, we show that the y-lines it induces are always accompanied by 
a continuum flux of cosmic rays. Hence constraints on cosmic-ray fluxes imply constraints on the intensity of 
y-lines and vice-versa. A comparison with up to date observational bounds shows the possibilities to observe or 
exclude cosmic rays associated to y-line emission, so that one could better determine the properties of the DM 
particle, possibly discriminating between some of the operators. 



I. INTRODUCTION 

One of the best "smoking-gun" signals for establishing the 
existence of an annihilating or decaying dark matter (DM) 
particle is the possible observation of a cosmic y-ray line HI. 
Forthcoming satellites |2] and airCherenkov telescopes |I3}{5], 
but also current instruments like the Fermi large area telescope 
IE 13 and the HESS instrument [8' "9], will allow to probe 
this possibility with further sensitivity. From DM particle an- 
nihilations the amount of monochromatic y-rays produced is 
expected to be limited (even if in some cases it can saturate 
the present observational bounds) because it is in general loop 
suppressed with respect to the total annihilation cross section, 
that is generically constrained by the DM relic density. For 
a decay, on the contrary, the amount of monochromatic rays 
emitted could a priori be much larger. Even if in many scenar- 
ios the DM particle is not expected to decay at all, there ex- 
ists a well-motivated theoretical framework where DM would 
naturally decay with a lifetime larger than the age of the Uni- 
verse: if its stability is due to an accidental low energy sym- 
metry that has no reason to be respected by any ultraviolet 
(UV) theory, just as expected for the proton in the standard 
model (SM). In this case the effect of the UV physics causing 
the decay is suppressed by powers of the UV scale. At low 
energy, this can be parametrized in full generality by writing 
down the most general effective theory respecting the low en- 
ergy symmetries of the model. For a DM particle mass around 
the electroweak scale (as supported by thermal scenarios), it 
turns out that a decay suppressed by 4 powers of the grand 
unified theory (GUT) scale, i.e. induced by dimension 6 op- 
erators, leads to a lifetime which can give fluxes of cosmic 
rays (CRs) of the order of current observational sensitivities 
ifTol [TTl . This is a too nice opportunity to probe the GUT 
scale to not study it in more detail. For an explicit example of 
a setup with an accidental symmetry, which can lead to DM 
decay into y-lines through dimension 6 operators, see refer- 
ence 1 1 r|. 

The use of an effective theory picture for a DM decay is 
fully justified because, unlike an annihilation, it necessar- 
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ily requires a large UV scale (unless one would invoke ex- 
tremely tiny coefficients). As explained below, this model- 
independent gauge invariant effective operator approach im- 
plies that when a DM particle decays, it does so to several 
final states. It is therefore different from the usual model inde- 
pendent approach of considering indirect detection constraints 
on separated single decay final states, as e.g. Refs. |12|. The 
amount of CRs that these several final states induce is also in 
general larger than those from electroweak corrections on a 
y-line final state, as considered in Ref. |131. By allowing a 
systematic determination of the y-ray line emission possibil- 
ities, and by properly taking care of the low energy symme- 
tries of the theory, the low energy effective operator language 
is the appropriate one to approach a series of important re- 
lated questions: What are the theoretical expectations to see a 
y-line from a DM decay (including, what effective structure a 
UV theory must induce to do so)'|j What are the possibilities 
to see a y-line, given the CR constraints, and conversely to see 
CRs associated to the observation of a y-line? And from there, 
what are the possibilities to discriminate among the various ef- 
fective radiative operators, hence to have information on the 
DM particle whose decay would have produced the y-line, and 
maybe on the UV physics associated? 



II. FULL LIST OF RADIATIVE OPERATORS 

To write down the most general effective theory of DM de- 
cays is a quite tremendous task, especially if one allows for 
new particles into which the DM particle could decay. In sce- 
narios of accidentally stable DM, such particles can be ex- 
pected, see e.g. Ref. |TI]. However, for the production of 
monochromatic lines it turns out that this can be done, i.e. the 
number of operators (up to dimension 6) is quite limited. 

This is due to the stringent criteria an operator must fulfill 
in this case: 



For this question it is important to keep in mind that the radiative operators 
we will list below could be induced either directly from the UV physics, or 
from low energy loop correction to other UV induced DM decay operators. 
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a) First, obviously the operator must contain the DM field. 

b) Secondly, in order for the operator to lead to a two-body 
decay with a photon, it must not contain too many fields 
(except eventually scalars that could be replaced by their 
vacuum expectation values (vev)) and it must contain ei- 
ther a hypercharge Fy or a SU{2)l Fl" field strength. If 
the DM particle is neutral, as we will assume here, the 
photon cannot come from a covariant derivative because 
in the two body decays all fields are necessarily neutral. 
One could eventually have a photon emitted from an op- 
erator that contains a new U{\) gauge field that kinemati- 
cally mix with the hypercharge gauge boson, L 3 zFy'^F^^ 
(rendering the various neutral particles to be effectively 
milli-charged). We will not consider this possibility here. 
It gives a range of bounds on the emission of y-ray lines 
which is similar to the one obtained below (see Ref. llT4l 
for details). 

c) Thirdly, some operators can be related to other ones 
through various relations, equations of motions, shift of a 
derivative or use of the fact that the commutator of two co- 
variant derivatives D^, gives F^i^. However one must be 
careful in using these relations. The criteria we apply here 
is that through these relations an operator can be dropped 
from the list only if in this way there is a one-to-one corre- 
spondence between this operator and another one already 
in the list. Otherwise both operators must be kept because 
in general they give different ratios of y-line to CRs. 

Applying the criteria above there are only three possible 
general dimension 5 structures, one for each of the three types 
of DM particle we consider here, scalar, fermion or vector: 
^dmF^vF^'\ ^DM<^f.v^F"\ Fl^fF^'^ respectively. Specify- 
ing what is the F^^ field strength one obtains only 9 possible 
operators, 5 for a scalar candidate, and 2 for a fermion or vec- 
tor candidate 
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where (^dm/'^dm denotes the multiplet whose neutral compo- 
nent ^^qm/^dm 1^ particle. By "(«,}')" we specify 
what must be the size n of the SU {2)i multiple ts and their hy- 
percharge Y . F'^^ stands for a new possible low energy gauge 
field and the vector DM operator f£J^ stands for an abelian or 



non-abelian DM field strength (in practice it will not be nec- 
essary to make this distinction in the following). \(/ and (|) are 
meant to be either SM fields when allowed by gauge invari- 
ance or new low energy fields. The symbols A — Exj„.v stand 
for a classification of the operators' possible astrophysical sig- 
nals, and will be explained in Sec.|IV| 

As for the dimension 6 operators the number of possibilities 
is also remarkably limited. Two general structures are singled 
out for the scalar case and three for the fermion and vector 
cases, leading to 7 scalar operators 
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to 6 fermion operators 
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to 5 vector operators 
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By Dj] we mean a covariant derivative that contains the DM 
vector field0 

Along the operator criteria defined above, note that one 
must still add to the list three types of operators, whose struc- 
tures are somehow more involved. 

First, a few operators that contain two covariant derivatives 
on a same field. There are two structures of this kind for 



^ In Ref. fl 11 an explicit example can be found of an accidental symmetry 
setup leading to the operators of Eqs. \24) and Note also that in 

Ref. 115 1 there are examples of heavy scalar and heavy vector setups whose 
exchange induces dimension 6 four-fermion interactions that at one loop 
induce a \\Idm — > 7V decay. The effective amplitude for this process is the 
same as the ones that the dimension 5 operators of Eq. |6} or Eq. JTJ give. 
This exemplifies the fact, to keep in mind, that dimension 5 operators for a 
decay can naturally have a "dimension 6 suppression" of the lifetime. 
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the scalar case, DfjD^i^oM^F^'^ , (^DMDijDy(^F'^ and one for 
the vector case, D^*'D^*'(t)(j)'F''^ (with F>'^ the hypercharge 
or SU{2)i field strength). These operators are reducible to a 
combination of operators where the two covariant derivatives 
are replaced by field strengths. 

Secondly, there are four fermionic operators with a 

Vi'MO;uv^V^y'L and ^DMPf^fiv^Fy^L- Through the equation 
of motion of the spinors, these operators are reducible to the 
same operator with^ replaced by a mass or a scalar field or a 
combination of several of these operators. 

Thirdly, there are also operators with dual fields. The full 
list of operators with dual fields can be obtained by replacing 
one field strength in the operators of the list above by its dual 
one. Note nevertheless in this case that Eqs. (|6])-(|7]i and ( 17 i- 
( fTSj l do not lead to any new operator (since they give back the 
same operator with the non-dual field) and Eqs. ([T5]l-( 16 1 do 
not give any y-line. Also Eqs. ([T9]l-([22)i give operators that are 
reducible to a combination of the same operator without dual 
and of the operator of Eq. (|6]l or (|7]|. 

In the following we will consider the phenomenology of in- 
dividual operators, as usually done in effective theories. Few 
comments on the possible effects of considering linear com- 
bination of several operators can be found in Sec. |VII| To 
discriminate operators, or sets of operators, it is not sufficient 
to observe a y-line at a given energy and with a given intensity, 
because these observables can be accounted for by any opera- 
tor, by adjusting the DM mass and the 1/A or 1 /A" coefficient 
in front of the dimension 5 or 6 operators. Additional informa- 
tion is needed, either on the y-line spectra they give (Sec III 
or on the associated continuum spectrum of CRs (Sec. IV VI 



Note that in the following we will not consider explicitly the 
phenomenology of operators above that are combinations of 
the operators of Eqs. (Tp7 i, i.e. of operators with 2 covariant 
derivatives on the same field, of operators with alp and of the 
operators of Eqs. ([T9]l-([22]i with dual fields. As for all other 
dual field operators, they give the same phenomenology than 
their non-dual partner. 



III. SINGLE OR MULTIPLE y-LINES? 

Let us start by discussing a simple way to, in principle, 
distinguish between some operators in the above list: from 
the different number of y-lines they can give and the relative 
strength and energy of these y-lines. Fig.[T]shows all types of 
two-body final states including y-lines that Eqs. ([T)-([27]i give. 
The operators giving more than one y-line can be separated 
into two categories, depending on if they include one or two 
SM field strengths. 

On the one hand all the operators with two SM field 
strengths in the operator give peaks associated to the yy 
and yZ channels. This could be the case for six operators, 
Eqs. (1-3, 10-12). The relative strength of the two y-lines are 

, cos^9iy sin^29H' sin^Giy 

2sm^eiv 2cos^29w 2cos^9w 




-,/z 



7, Z, Z\ 0, 0' 



FIG. 1. The y-line decays the operators can give and their associated 
non radiative decays. 



F,^^ and /i^v^f^ re- 



for the operators including Fy, 
spectively, with F = (wdm ^ "^z)^ ^^'^ Weinberg 
angle. If itidm is above but not too far above the niz scale, 
these three types of operators could thus be distinguished 
from each other because both peak strengths could be sepa- 
rately resolved at energies equal to mom/'^ and moyv//2(l — 
m|/m^^)). In practice only one y-line will appear when the 
detector resolution and sensitivity are not good enough to sep- 
arate too close monochromatic lines. 

On the other hand there exists also 6 operators with only 
one SM field strength which could also give two y-lines: the 
scalar DM ones in Eqs. ( 15 1-( 16 1 which allow a yZ and yZ' 



decay (with a Z' from the covariant derivatives), the vector 
DM of Eqs. (|24]i-((25]l and Eqs. ((26])-(|27]) which if (|) ^ (|)' 
allow both a y(|) and a y(|)' channel. If the mass of the Z', 
(|) or (|)' differs sizably from mz these operators would typ- 
ically give different energy gaps between the y-lines com- 
pared to the operators with two SM field strengths. In par- 
ticular for these operators, if these masses are not many times 
smaller than moM, and have sufficiently small decay widths, 
one could have an observable multiple peak structure also for 
moM > mz; at moM/^il - m^/m^j^j) ~ mdm/'^ energy and 
at ot/)m/2(1 — m^i ^ ^//'Wom) energies. That would single out 
these 6 operators. Information on the height and energy of 
both peaks provides further informations on the masses and 



couplings involved in the effective operators. For Eqs. ( 15 1 
and ([T6| one has 
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with g the SU{2)i coupling and gy'^Y' are the extra gauge 
group coupling and charge of the scalar field under this gauge 



group. For Eqs. (24 1-(27 1 one has 
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with ((])) and ((|)') the (|), (|)' vevs. None of these 6 operators give 
any yy final state. 
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The remaining 15 operators, i.e. all the fermion DM par- 
ticles' operators Eqs. ([6]- [9|17f|22|, and the scalar and vector 
DM particles' operators in (4 - [5|13 - 14 23 i, give only a single 
y-ray line, and would thus be excluded as an explanation to the 
observation of multiple monochromatic y-lines. Of course, a 
DM particle can always give multiple lines by appearing in 
several operators involving fields with different masses. 



IV. CONTINUUM RAYS ASSOCIATED TO A y LINE 

In the previous section we have seen how multiple y lines 
might allow to single out some operators. We will now fo- 
cus on a feature shared by all operators in our list, which will 
allow us to characterize them further: all allowed operators 
that give a decay into a y-ray line will unavoidably also pro- 
duce, with different amounts, decays into particles producing 
CRs. This is a consequence of the fact that gauge symmetries 
dictate the possible effective operator structures. It is particu- 
larly clear that, if kinematically allowed, a Z boson is always 
among the final states since the hypercharge and SU{2)i neu- 
tral gauge eigenstates are made of a mixture of a photon and 
a Z boson|3 

This multiple final state property is the crucial one why, by 
merely studying the corresponding continuum of CRs at lower 
energies, we can put upper bounds on the possible strength of 
the y-line signal for any given operator. In the following we 
will focus our analysis on the large nioM > 200 GeV mass 
region, where the production of CR from a given operator is 
unavoidable, i.e. nioM is sufficiently large so that channels 
with up to two SM gauge bosons are kinematically fully open. 
We will also make the assumption that the mass of any new 
beyond SM final state particle is negligible with respect to the 
DM particle, so that if the operator leads to several peaks, they 
all show up at the same Ey ~ mj)M/2 and their contributions 
are summed|3 

Inspection of the full list of operators shows that there is 
no single operator which can give a ratio of monochromatic 
photons to Z production larger than (what we call prediction 
A) 
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where CR ~ p,yc,e^,... stands for any CR which could be 
produced out of the Z, with nc«/z the number of CRs pro- 
duced per Z boson at a given energy. In contrast to the 



^ This also holds for a DM annihilation scenario when the 7-line(s) is induced 
by effective operators [16t'18'| or induced at one loop in explicit models 
11911211 . This property is also characteristic of specific gravitino decay 
scenarios ( 22 1 . 

If on the contrary the mass of this new state is close to moM, some chan- 
nels could be more phase space suppressed than others, modifying the con- 
straints. Note also that, if one decreases mou towards the mz threshold, 
the Ry/CR r^tio could go to zero or to infinity, depending on whether the 
operator induces a decay to yy. 



monochromatic y-line, the y^ stands for the continuum of 
lower energy photons produced by hadronization and elec- 
troweak correction ||23]| connected to the production of Z 
bosons|^ Given the fact that CR fluxes are observationally 
bounded from above, this ratio sets the absolute upper bound 
on the intensity of y-ray lines that a single operator can give, 
to be determined below. 

It is out of the scope of this work to list all the exact predic- 
tions the various operators give for the monochromatic pho- 
tons to CR emission (and the corresponding upper bound they 
give, see Ref. 1(141 ) because a single operator may give differ- 
ent ratios depending on the quantum numbers of the particles 
it contains, and the nature of the (|) or \\t particles they might in- 
volve. However the A — Ex,m.v labels we have put beside each 
operator in Eqs. ([T]i-([27]i allow already to see rather well what 
is going on for each of them. Varying between all possible 
multiplets in the operators, each operator turns out to give a 
maximum ratio which is given either by A above or by one of 
the four following relations 
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(32) 
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with — 1/4, 1 for D and E, respectively. Signal ratio A is 
characteristic of operators containing only F^^, B is typical 
of those containing both F^'^ and F^^, and signal ratios C to 
E correspond to operators including only F^^ . These bounds 
are obtained assuming that the (|) and \\t particles which may 
show up in the two-body final state do not produce any CRs 
(as could be the case for a beyond SM particle). 

From these maximal ratios, each operator may give a Ry/cR 
ratio that can be smaller in three possible ways, if the extra 
(|) or i|/ particle does decay to CRs (in these cases we put a 
"x" index label), if the operator gives different ratios for dif- 
ferent multiplets ("m" label), or if the Ry/cR depends on the 
(unknown) vev of a scalar field ("v" label). There are also 
operators which give a totally fixed ratio. In these cases we 
do not attach any index label. This is the case for two A type 
operators and for all B cases. 

The "x" and "m" variations of Ry/cR that can arise for an 
operator are relatively limited. Given the fact that in each case 
there is already production of CRs from the Z (and sometimes 
W) boson, if we assume the extra particles "x" give an amount 
of CR similar to what a Z or a W boson produces, these par- 
ticles cannot modify largely the ratio. Typically, the change 
is not more than a factor of a few or at most, in rare cases of 
type A predictions, around one order of magnitude. Similarly, 
among the different multiplet configurations that do lead to a 



^ Final state radiation and internal Bremsstrahlung processes could mimic a 
line-like signal. We do not account for these processes in the numerator of 



. ratio. 
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y-line, a given operator can have the RyjcR ratio varying by a 
factor of a few up to at most a factor 5 — 6 (except for Eq. ( fTS] ! 
which allows a factor ^ 10 variation because it involves two 
fields on top of the DM and gauge boson fields). For exam- 
ple many operators of class C with an index "m" may easily 
give, for some multiplets, bosons in the final state, in that 
case the signal ratios can give predictions equal or close to D 
or E. Therefore from these effects it is not possible to go far 
below the maximum signal ratio, A-E, that each operator are 
assigned. Only by varying the value of the vev of a beyond 
SM scalar field can we go far below the operators maximal 
Ry/CR ratio value, which is only possible for the two "v" la- 
belled operators in Eqs. 



15 1 



16 1. The latter possibility stems 
from the fact that the Z + (|) production is proportional to m| 
whereas the y-line productions are instead proportional to the 
vev of (|), and thus disappear when the vev goes to zeroj^ 

As examples of "v" specific cases which necessarily give 
ratios much smaller than the A-E predictions, let us consider 



iSf being a triplet of hypercharge 2 in the operators of Eq. ([T5 1 
and 



16 1 



We call these predictions Fy and F^ respectively. 
For these setups the y-line production is proportional to v|, 

whereas the Z(|) production is proportional to Since the 
vev of a triplet is constrained to be small by electroweak data 
(|v| < 4 GeV), the R^/cr is suppressed by at least a factor 



vl/m| ^ 10 ^. As aresult, for example the Fi case of Eq. (|16 1 



lies a factor ^ 1000 below the absolute one of case A forme 
RylCR ratio (and with a y-line intensity so weak that basically 
it cannot be distinguished from the associated continuum of y 
rays, see below). This shows that operators that give "v" la- 
belled predictions crucially depend on the electroweak quan- 
tum numbers of the particles involved and for some multiplets 
predict very weak lines. As stated above the same operators 



of Eq. ( 15 1 and Eq 



16 1 can saturate predictions A and C (for 
instance for isf a SM singlet, whose vev can be large and break 
an extra gauge symmetry, so that the DM particle can decay 
to the associated gauge boson(s) and a y). 



V. GENERAL RESULTS 

Once the ratio between the predicted line signal and the 
lower energy energy continuum of CR is known, observa- 
tional constraints on CRs put indirect limits on the possible 
strength of any associated y-line signals. In Fig.|2]we show for 
the predictions A-E what are the corresponding upper bounds 
on their decay rates Fy into monochromatic y-rays, imposing 
that their continuum-part of CRs are not overproduced. The 
strongest indirect bounds we get on the y-ray lines are from the 
R4MELA anti-proton data f24] (dashed lines) or the Fermi- 
LAT isotropic y-ray flux measurement 125] |26ll (solid lines), 
depending on the value of nioM- 

The limits are obtained by taking = ^,(data,- — 
model,)^/a? = 9. Here "model," refers to the sum of the DM 
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FIG. 2. Upper bounds on the decay rate into monochromatic pho- 
tons for the predicted ratios A-F. Dashed curves represent the indi- 
rect constraints imposed by antiproton data from PAMELA, whereas 
solid curves are derived from the isotropic diffuse y-ray measurement 
by the Fermi-LAT. The grey shaded area represents the direct obser- 
vational exclusion region from Fermi and HESS data. 



signal and our background flux model, and the standard devia- 
tions o, are taken directly from |24| and [26l. The anti-proton 
flux can be well described by conventional astrophysics, and 
above 10 GeV (where solar modulation has a small effect) the 
canonical setup has relatively small uncertainties in its pre- 
diction. We therefore use only data above 10 GeV, and take 
our anti-proton background-flux model to exactly match the 
observation in each energy bin {j} = by construction). The 
origin of the isotropic diffuse y-ray background is still quite 
uncertain and we thus take that background flux to be a com- 
pletely free (positive definite) function of the y-ray energy]^ 
For the DM signal we take a NFW DM distribution (as in 
f23|) and include the extra Galactic contribution using the op- 
tical depth from Ref. [,27 J . The Galactic signal is not fully 
isotropic, and we make the reasonable choice to consider the 
DM signal from the Galactic poles {\b\ ^ 90°). For the anti- 
proton flux we choose the MIN propagation model (defined 
in Il23ll28l ) to always stay as conservative as possible in our 
below statements. For the MED and MAX propagation mod- 
els (defined in ll23l |28]| ) the various A-F curves in this figure 
would be approximately suppressed by a factor of 4 and 10, 
respectively. As Fig.|2]shows, between the A and E cases the 
y-line constraint varies by as much as a factor ^ 50. 

In Fig. |2] these indirect upper bounds on Fy can be com- 



^ In practice, this means we only consider the y flux from DM in bins where 
* Note that for Eqs. j24} - j27| there is also a vev dependence but it is much it overshoots the observed flux when calculating Ax~ (as y} = Ohy con- 

milder, struction in the remaining bins). 
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pared to the current direct exclusion bounds on y-ray line 
signals (grey shaded areas). Up to Ey — 200 GeV we used 
the limits from the y-line searches by the Fermi-LAT collab- 
oration |6|. From 200 GeV to 500 GeV we used the same 
— 9 procedure as described above, applied to a y-line on 
top of the isotropic diffuse y-ray measurement ll26l . Above 



Ey = 500 GeV and up to 25 TeV the line signal constraints 
from HESS [8 1 are translated to limits on decaying DM. When 
using the HESS collaboration y-line flux limits, we derived 
constraints by assuming a NEW DM profile (the limits are 
similar for an Einasto profile). We studied both their isotropic 
sky and central Galactic halo region data, but retain only the 
strongest constraints. 

From Fig.|2]one can directly see what could be the implica- 
tions of a detection in the near future of a y-line, of a CR ex- 
cess or of both. To this end it should be kept in mind that if the 
operator considered carries a "x", "m" or "v" label, depend- 
ing on what are the exact fields considered in this operator, its 
corresponding bound can be suppressed by factors discussed 
above compared to the maximal ratios Ry/cR shown in Fig.|2] 

First of all, as Fig.|2]shows, the A-E ratios turn out to have 
y-ray line upper bounds from their CR constraints which are 
weaker than the present direct observational limits on y-lines 
(in a mild way for E though). This has the important impli- 
cation that, with today's CR constraints, any of the operators 
could accommodate a possible near future situation where a 
y-line would be observed without any observation of a CR 
excess. Only special cases could be excluded along such a 
scenario; basically the "v" labelled operators for given config- 
urations such as the Fy and Fl upper bounds shown in Fig. [2] 
(and marginally also some special "x" or "m" suppressed C-E 
cases). Conversely, a broader energy continuum of photons, 
just below current constraints, could not be due to operators 
with the A to D predictions. Similarly, the observation of an 
anti-proton excess (using in this case the "MAX" propagation 
parameters to stay conservative) could not be due to operators 
with predictions A to B. Obviously, an anti-proton (gamma 
continuum) excess can neither be attributed to an operator's 
prediction if the DM particle mass is in a range above (below) 
^ 1 TeV where the y-continuum (anti-proton) constraints are 
stronger 

This pattern could change in the future depending on with 
which precision a CR excess and/or a y line would be observed 
or bounded. With the upcoming AMS-02 data, and the po- 
tential to restrict the charged propagation model to be closer 
to the MED or MAX propagation models, the CR continuum 
constraint could potentially be improved by more than an or- 
der of magnitude for nioM below a few TeV 1291 . In this case 
the A-E lines would lower accordingly. Similarly, Gamma- 
400 (2 1 and/or the CTA telescope 1 3 -5 1 could improve the y- 
line sensitivities by one order of magnitude depending on the 
energy range and sky region (see fl5J for the case of decaying 
DM). 

For instance if a y line is observed, then the ratio of the ob- 
served y-line decay rate and the maximally allowed y-line de- 
cay rate, as given by the (updated) A-E predictions of Fig. [2] 
will tell us which operators can be excluded and which one 
could lead to an observation of an associated CR excess. If 



this ratio is larger than 1 the operator is excluded (and of 
course even more excluded if there are "x" or "m" suppres- 
sions, as this would increase the CR signal for an operator). 
If it is smaller than one, but not much smaller than one, then 
the CR excess associated to this y-line is expected to be de- 
tectable. 

Put in another way, taking into account these possible future 
experimental improvements, the following picture emerges. 
On the one hand, for the operators that involve the Fy"^ field 
strength, A and B cases, and if they have no "m" or "x" or "v" 
possibility of suppressions, it appears clearly improbable that 
in a near future we could see a CR excess associated to a y line. 
For A cases with index "m" and/or "x" this is not impossible, 
but difficult. For the A case with "v" label, it is possible to see 
both associated signals easily, even without invoking any "m" 
or "x" suppressions (as in the Fy example). On the other hand 
for the operators that have a field strengths, C-E cases, 
the discussion is the same but with all upper bounds shifted 
downward, so that it is easier to see both associated signals. 



VI. STANDARD MODEL FINAL STATE RESULTS 

We now focus our analysis on the special case where all 
particles beside the DM one would be SM fields. This turns 
out to be a possibility for all operators (with \\t replaced by the 
e, /J or z doublet and the scalar (|) fields replaced by the SM 
scalar doublet H = {H^\H^Y or its conjugate), with the ex- 
ception of Eqs. ([8]), (|9]l, ([13]), ([14]| and (|25]|. Note that there is 
no dimension 5 operator of this type for vector DM particles, 
and for operators of dimension up to 6 the scalar DM parti- 
cles can only have the quantum numbers (1/3/5,0), (2/4/6,1), 
(4/6,3) or (6,5). For fermion DM particles only the (1/3/5,0), 
(3/5,2) or (2/4,1), (4,3) multiplets are possible. The values of 
Ry/CR when the operators include only SM decay products are 
particularly interesting because they are almost totally fixed. 
Indeed, in these cases the Ry/cR ratios neither depend on any 
unknown final state that might or might not give CRs, nor on 
the unknown vev of new fields. For a given operator the RyjcR 
ratio may only depend on the way the SU{2)l indices are con- 
tracted. Here we will consider only the case where the DM 
particle has vanishing hypercharge {e.g. to avoid constraints 
on direct detection processes mediated by a Z boson) and the 
so called "inert doublet" candidate where the DM multiplet is 
a scalar doublet with the same hypercharge as the SM scalar 
doublet. 

Depending on the operator, the hyperchargeless candidates 
lead to predictions A-E or new predictions, which we call 



n F 



A/j and C/,. The predictions Z)g^ x and /if.^ x 



correspond to the cases where a e, /j or T charged lepton is 
produced in association to a W in the final state. They are 
obtained from the D and E, respectively, by simply replacing 
«c«/w± by «CR/w± +ncR/i^ in Eq. (|34]). The predictions A/, 



and C/,, corresponding to when a SM scalar boson is produced 
with the y and Z in the 2-body final state, are obtained by sim- 
ply adding a nc«//! term in the denominator of the prediction 
A and C, respectively. The resulting constraints are given in 
Fig.|3] and the dictionary between signal predictions and pos- 
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Signal R^icR 


Operator in Eq. 


Quantum no. {n,Y) 


A 


(1) 


feM^l-O) 




(17) 


¥dm=(1,0), (3,0) 


Ah 


(24),(26) 


Vdm = {1,0), HH = {1,0) 


B 


(2) 


<^DM=(lfil (3,0) 


C 


(18) 


Vdm=(3,0), %H = {1,0) 


Q, 


(25),(27) 


Vdm = {hO), HH = (3,0) 


D 


(18) 


WDM =(5,0) 


D, 


(18) 


Vdm =(3,0), xj/DM^ = (2/4, 0) 


E 


(18) 


WDM =(1,0) 


El 


(18) 


Wdm =(3,0), WiWdm = (2/4, 1) 


no-lines 


(3) 


<t)DM=(3,0) 




(18) 


Vdm =(3,0), V/^/ = (1/3,0) 



DM 



[GeV] 



TABLE I. All possible signal scenarios for a hyperchargeless 
DM candidates decaying into SM particles. The labels of the 
monochromatic-line to cosmic-ray signal ratios Rj/cr are explained 
in Eq. l[3TJ - ^ and in the text of Sec.jVl] 



FIG. 3. Same as for Fig.[2]but for SM final states. Note that only the 
Ee.fj.T and not the D^^^ x predictions have been plotted for readability 
sake, but their ratios are the same as the EXoD ratio that can be found 
inFig.|2] 



sible effective operators for these hyperchargeless multiplets 
is given in Table [l] 

For the inert scalar doublet case one can get A, B, E or 
new len gthy pred ictions llT4l from the "v" labelled operators 



161), which we call g(/'^ and G^' and depend 



I/IR 



in Eqs. |l5| - 

on whether the DM particle is made of the real (R), of the 
imaginary (I), or of both components (RJ) of the neutral com- 
ponent of the inert doublet (taking the SM scalar vev to be 
real by convention). Depending on whether the A-5 term in the 
tree-level scalar potential gives a positive, negative or vanish- 
ing mass splitting between the two neutral components (see 
e.g. EOl), these G'l"^ predict ions are different because, as for 
the A,5 term, these two operators break the Peccei-Quinn type 
symmetry between these two components. The signal to op- 
erator dictionary for this case is given in Table [n] 

Fig.[3]shows that, similarly to what one observed in Fig.|2] 
almost all cases with only SM final states have potential y-ray 
line signal strengths that are not restricted by their associated 
continuum of CRs. The only exception where actual direct 
limit on a y-line are weaker is for the Gl cases - mainly be- 
cause the p hoton-line emission for this "v" labelled operator 
in Eq. ( 16 1, is suppressed by a factor of = 1/4 with respect 
to the W^W^ emission. 

Note that for the cases with prompt charged leptons in the 
final state, leading to {D/E)e^^.z predictions shown in Fig. [5] 
we have checked that the y flux from inverse Compton inter- 
actions of these leptons can be neglected. As for the observa- 
tional positron and electron constraints ||32l - |34l . they turn out 
to not give any further bounds here. 

Note finally that in Table. 1 there is a case where the DM 
field is a quintuplet. Assuming that there is no beyond the 



Signal R^icR 


Operator in Eq. 


Quantum no. (n, Y) 


A 


10 






B 


11 






E 


12 




FlFl = (1/5,0), iSfDMFL = (2/4,-1) 


I/IR 

(jy 


(15 




DM^Im[*Vl/Re+Im[<J 




(16 




DM^Im[*V]/Re+Im[*V] 


no-lines 


12 
15 


,(16i 


FlFz, = (3,0) 
DM^Re[<^] 



TABLE II. Similarly to Table I, all possible signal scenarios for an 
inert scalar doublet = (2, 1) decaying into SM particles. 



SM fields other than this quintuplet, this is the quintuplet of 
the "minimal DM" setup 131]. The neutral component of this 
quintuplet is naturally stable enough because "accidentally" it 
cannot appear in any operator of dimension less than 6. Ta- 
ble. 1 shows that it can lead to a y-line through the operator of 



Eq. (18 1, leading to prediction Eg^^j, depending on which SM 
lepton doublet is considered in the operator, see Fig.|3] 



VII. A FEW FINAL COMMENTS 

In Figs. [2] and [3] we took into account only the full list of 
two-body decays the various operators give. This is expected 
to be a good approximation, except for the following case: 
if the two-body decays proceed through the vev of a scalar 
and if the DM mass is far above the value of this vev. In 
this case, if the three-body decay is kinematically allowed, 
the two to three-body decay-width ratio is expected to be of 
order ^ 487t^((|))^/OT|)^ . Finally note that a UV theory is of 
course not necessarily expected to induce only one relevant 
operator In presence of several operators the CR production 
could be easily boosted relatively to the y-line signal, so that 
the A to G y-line bounds in Figs. |2] and [3] would lower accord- 
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ingly. Alternatively, with several operators there exists also 
the possibility that the UV physics induces several effective 
operators that would produce the same final states, leading 
to a destructive or constructive interference for some of the 
channels. This issue is in many respects similar to the well- 
known one of determining the cut-off of possible new physics 
that would induce operators composed of only SM fields f 35l 
(i.e. lower bounds on UV cut-offs given in the literature are 
typically derived for single operators, and they could be de- 
creased if interferences between operators are allowed to be 
tuned). The variation in Ry/cR that such interferences could 
cause does highly depend on the operators considered and the 
multiplet considered in them. For many cases it is not possi- 
ble to increase sizably Ry/cR in this way. For the other cases it 
is possible by invoking one or several tunings of the operator 
coefficients (i.e. of couplings and quantum numbers in the UV 
theory). 

In a similar vein one has also to keep in mind that, through 
renormalization group effects, a single operator generated at 
a high energy scale can generate other ones at low scale (al- 
though in general with sub-leading coefficients). 

In conclusion only a few dim < 6 effective operator struc- 



tures can lead to a monochromatic y-line from DM decay, even 
if one accounts for possible new BSM particles in the decay 
products. Each of these operators leads to a y-line to CR flux- 
ratio which turns out to be bounded from above and in almost 
all cases also from below. This leads to clear possibilities of 
observing a continuum of CRs associated to a detected y-line. 
From the observation of a y-line and/or of a CR fluxes that 
these operators could induce, and from the observation/non- 
observation of multiple y-lines (with characteristic energy and 
height), there exists a clear potential to constrain/discriminate 
various sets of operators. 
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